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Abstract 



"The energy deposition and redistribution in hot Jupiter atmospheres is not well understood currently, but is a major factor for their 
i evolution and survival. We present a time dependent radiative transfer model for the atmosphere of WASP- 18b which is a massive 
C> (10 Mj up ) hot Jupiter (r eq ~ 2400 K) exoplanet orbiting an F6V star with an orbital period of only 0.94 days. 

Our model includes a simplified parametrisation of the day-to-night energy redistribution by a modulation of the stellar heat- 
er ing mimicking a solid body rotation of the atmosphere. We present the cases with either no rotation at all with respect to the 
<D synchronously rotating reference frame or a fast differential rotation. 

The results of the model are compared to previous observations of secondary eclipses of iNvmever et al.1 (1201 II) with the Spitzer 



t — Space Telescope. Their observed planetary flux suggests that the efficiency of heat distribution from the day-side to the night-side 
of the planet is extremely inefficient. Our results are consistent with the fact that such large day-side fluxes can be obtained only if 
there is no rotation of the atmosphere. 

p I , Additionally, we infer light curves of the planet for a full orbit in the two Warm Spitzer bandpassses for the two cases of rotation 
fT"] and discuss the observational differences. 

Keywords: Atmospheres, structure, Extrasolar planets, Radiative transfer 
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1. Introduction 



WASP- 18b jHellieretalll2009l) is a massive planet orbiting 
its F6V parent star. The high mass and short period of WASP- 
18b make it a useful test-case for explori ng the tidal interac- 
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tion between pla n ets an d their host stars (IBrown et al 
iMatsumura et al. . 2010h . Accurate parameters for WASP 



18b from high quality lightcurves have been presented by 
ISouthworth et al. ( 2009 ). The observed Rossiter-McLaughlin 
effect for WASP- 18 suggests that the orbi tal axis of the plane t 
is aligned with the rotatio n axis of the star Triaud et al. ( 2010h . 



Nvmever et al 



(1201 lh studied two secondary eclipses of 
WASP- 18b obtained with the fourlRAC channels of the Spitzer 
Space Telescope. As a consequence of the high mass of the 
planet, the signal we measure originates from deeper in the at- 
mosphere (higher pressure) than for a lower mass planet. These 
author infer indeed a contribution function peaking at 10 bars 
for 3.6 pm in a case without thermal inversion. For a lower 
mass planet such as HD 189733b (M p ~ 1.14M Jup ), the same 
bandpass sounds an atmospheric pressure of around 0.3 bar 



dKmitson et alll2009j) 



The analysis of INvmever et al.l (1201 lh suggests that the 



planet has a very low albedo as well as an inefficient heat re- 
distribution. They present ed a parametrised atmospheric m odel 
following the method of iMadhusudhan & Seagerl (|2009) that 
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marginally favors a thermal structure presenting an inversion, 
as expected for plane t s receiving such high level of irradiation 



(Bur rows et all l2006t iFortnev et all 120081) . even though a non 



inverted model fits the observations within the error bars. 

This planet is among the ones we should be able to derive 
strong constraints on its heat redistribution or circulation pat- 
tern. Indeed, in principle, the Spitzer Space Telescope could 
be used to construct a day t o night map of heat r edistribution 
following the method from iKnutson et alj d2007l) . These au- 
thors observed HD 189733b during a full orbit in order to get 
the phase resolved brightness temperature distribution. This 
kind of study for this massive planet would enable us to test 
atmospheric models under a different regime by constraining a 
complementary pressure range. 

In Section|2] we describe the radiative transfer model used in 
this work and the method used to parametrize the heat redistri- 
bution. In Section[3] we present and discuss the results includ- 
ing the thermal structure and spectra as well as phase curves of 
the planet. 



2. Atmospheric model 

2.1. Radiative transfer 

We use line-by-line radiativ e transfer calcul ations from 0.3 
to 25 pm, described in detail inllro et al.1 (120051). 



In their modelling of WASP- 18b . INvmever et al. (120111) 
applied the method of Madhusudhan & Seagerl ( 2009 ) where 



10 7 models are computed with six parameters describing 
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the three layers of the P-T profile and four parameters control- 
ling the departure from chemical equilibrium of the four major 
species (H2O, CO, CH4 and CO2). A^ 2 minimization is per- 
formed between the numerous models against the four observa- 
tional data points. Their goal is to eliminate ranges of parame- 
ters that are too far from the observational constraints. 

On the contrary, in our approach we calculate the chemical 
abundances by assuming that the atmosphere is in thermochem- 
ical equilibrium with a solar abundance of the elements. In a 
more complete model, NLTE processes should be taken into ac- 
count. The model has 56 pressure levels from about 3.5 kbars 
to 10~ 3 mbar. By comparing the radiative lapse rate obtained 
with the adiabatic one, we adjust for the convective zone. 

As opacity sources, we include the main molecular con- 
stituents: H 2 0, CO, CH 4 , Na, K and TiO; Collision Induced 
Absorption by H2 and He; Rayleigh diffusion and H bound- 
free and H 2 free-free. We did not include CO2 but in thermo- 
chemical equilibrium, its abundance is low and we expect that it 
does not impact on the radiative budget. Moreover, the current 
model does also not account for clouds. 

We can calculate the evolution of the temperature pro- 
file in an atmosphere under hydrostatic equilibrium by using 
a time-marching algorithm, that allows us to include time- 
varying conditions. By using this approach, different from 
Madhusudhan & Seageri (120091) . we are able to introduce a 
parametrized prescription of a rotation of the atmosphere. This 
can be related to an estimate of the wind speed and, on the first 
order, the efficiency of the heat redistribution from day side to 
night side. 

2.2. Time-dependent calculations 

After convergence to a static solution corresponding to con- 
ditions where the planet is under globally-averaged insolation, 
we introduced a variation with time to the heating due to the 
star insolation. We set as outer boundary condition a time- 
dependent incoming stellar flux (more details can be found in 



dependent incoming s 
Iro&Demingll2010l) : 



F™(t,A,0) = 



4*-) cos [0(f)] F ( i(A) 
a{t)j 



(1) 



At each time-step, we calculate the star radius the star-to- 
planet distance a{t), 6 the planetary longitude of the point that 
is substellar at t = and F® the flux emerging from the star 
for each wavelength A. As a consequence, each parcel of atmo- 
sphere receive a stellar flux as if the global atmosphere were in 
rotation at a velocity set up by 6(t). The planet itself is assumed 
to be in synchronous rotation. 

We performed calculations in two cases: one where the in- 
coming flux depends only on the longitude 9, constant with time 
in that case (flux maximum at the substellar point, nill for the 
night side); and one where we introduce a solid body rotation 
(6 varies with time) corresponding to a velocity of about 6.5 



km s 1 at 1 bar. In Showman & Polvanil ( 201 1 ). superrotat 



ing jets at the equator with speed of the order of the km s _1 
can be maintained by angular momentum pumping from re- 
gions outside the jet to within the jet ('up-gradient eddy flux' 



proposed by 2008). Several eddy transport or planetary 
scale waves can be responsible for such transfer. A possibil- 
ity proposed by these authors is Rossby waves. It should be 
noted that ot her mechanisms c an generate this equatorial super 
rotation (see ICho et all 120081) . Additionally, in the radiative - 



hydrodynamical simulations by iDobbs-Dixon & Agoll (120121) . 
zonal winds develop with velocity from ~ 2 km s at 0. lbar to 
~ 4 at 10 3 bars, in the case of HD189733b. We then choose 
6.5 km s _1 as an exploratory value in order to asses an upper 
limit of the wind velocity. 

2.3. Input parameters 

The planet, star an d orbit parameters are taken from 
Southworth et al. I d2009h . They are summarized in Table Q] 



Table 1: Planetary, ste llar and orbital parameters used in the model (from 
Southworth et al.. 2009). 



Parameter 



Value 



Orbit 



0.02 AU 





P 


0.94 d 




e 


0.03 


Planet 








R P 


1.165 flj U p 




M p 


10.4 M Jup 




T 


2400 K 




Wday) 


3100 K 


Star 








Type 


F6V 




T e f{ 


6400 K 






1.23 R@ 



The equilibrium temperature is defined as: T eq = (1 - A)'^ 4 x 
/ 1/4 x x (R+la) 111 where A is the planet albedo; f the heat 
redistribution factor of the planet's atmosphere; r* the star's ef- 
fective temperature; the stellar radius and a the semi-major 
axis of the orbit. Using the parameter of the system from Ta- 
ble [1] this gives in the case of A = 0: r eq = 2400 K consid- 
ering a heating averaged over the whole planet (/ = 1/4) and 
2900 K considering an average over the day-side hemisphere 
(/ = 1/2). This makes WASP-18b a very hot Jupiter. For such 
highly irradi ated planets, we expect a temperature inver sion on 



2006; Fort nev et al.U2008l) 



the day- side dBurrows et al. 

The stellar flux has been computed by scaling an F- 
type star synthet i c tem plate from Castelli and Kurucz 
dCastelli & Kuruczl l2003lr l for an effective temperature of 
6400K. 

3. Results 

The temperature structure of the planet is illustrated by equa- 
torial cuts at the equator in Figure Q] The effect of the rotation 



2 See |f tp : //ftp .stsci ■ edu/cdbs/gr id/ck04models 



is obvious: the more efficient heat redistribution in the longi- 
tudinal direction leads to a more mixed atmosphere where the 
temperature contrasts are dimmer. The night-side gets warmer 
and the night-side cooler. In the no-rotation case, day-to-night 
temperature contrasts go from 2,000 K to 3,000 K depending 
on the pressure level. With a fast rotation, temperature contrasts 
drop to ~ 800K for pressures lower than 1 bar. 
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at ~ 2600 K and even the evening limb (90°) is hotter than the 
substellar point. These considerations are particularly interest- 
ing since it is the expected level that Spitzer observations are 
probing. 
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Figure 2: Thermal profiles for selected longitudes for the two cases considered. 
Upper panel: no rotation, lower panel: fast rotation. 



Figure 1 : Equatorial cuts of the planet thermal structure for the two cases con- 
sidered. Upper panel: no rotation, lower panel: fast (~ 6.5km s at 1 bar) 
rotation. The colors represents the temperature in Kelvin. The black line shows 
the direction of the star (substellar point). The pressure levels are shown in 
logarithmic scale with 10~ 3 mbar at the top of the atmosphere and one unit rep- 
resent one order of magnitude. For quick review, three white circles are drawn 
at respectively 1 mbar, 1 bar and 10 3 bar from exterior to center. 



Figure|2]shows thermal profiles at selected long i tudes. These 
can be compared to Figure 13 of iNvmever et al. I fcOllh . The 
deep (> 1 b ar) temperatures we cal culate are lower than the one 
obtained bv lNvmever et al.1 (1201 ll) . This is due to the fact that 
in our approach, we initially have the temperature averaged for 
the whole planet and not the day-side. Consequently, the net 
heat flux reaching the deep levels of their ID model is substan- 
tially higher than for the current work. In our model we have 
an inversion in the 1-0.1 bar pressure range. In the no-rotation 
case, the temperature is maximum at the substellar point peak- 
ing at 2900 K. When rotation is included, the thermal structure 
is more complex. An inversion is still generated on the day- 
side but the hottest profile is for a longitude of 45° with a peak 



Figure |3]represents the planet-to-star flux ratio as a function 
of wavelength for several positions - at the substellar point, 
at an intermediate longitude (representative of an average over 
day-side) at the night side and also for the planet in a planet- 
averaged stellar heating condition. 

Comparing the no-rotation with rotation cases, we can see the 
effect of decreased contrast (with rotation) between night and 
day. The cooler day-side exhibits indeed a lower flux whereas 
the warmer night side has an increased flux compared to the 
no-rotation case. 

The fast rotating case experiences flu xes too low on the day 
side to reproduce the values observed bv lNvmever et al.l(1201 ll) . 
whereas they are very well fitted in the non rotating case. This is 
consistent with these authors conclusion that WASP- 18b should 
have very little redistribution of heat from the day side to the 
night side. 

A trend seems to emerge: new observational constraints con- 
firm that very highly irradiated planets have very poor day 
to night energy redistribution, as another example is given 
by WASP- 12b whose phase curve has been measured by 
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Cowan et al.l (1201 lh . As stated by these authors, several fac- 
tors contribute to an ineffective redistribution for very hot plan- 
ets: magnet i c drag , very efficient at such large temperatures 



( Pernaetail 201 Op. as empha sized by the GCM simulations 



by iRauscher & Menoul (120 12l) where magnetic drag substan- 
tially slow down winds in HD209458b for a planetary mag- 
netic field of at least 3G (even if thermal ionization could be 
insufficient in nominal c ases, substantial ioniza tion by stellar 
radiation is possible, see Koskinen et al. , 2010t) : the radiative 
timescale decreasing with temperature a lot faster than the ad- 
vective timescale; and by the peaks of thermal emission and 
stellar heating being close, a weak greenhouse effect. We re- 
fer to the above articles for a more detailed study. One can 
also mention that h eat redistribution can h ave non-trivial mech- 
anisms by waves ( Watkins & Choi 2010l) and by instabilities 
Polichtchouk & Choi d20l3) 
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Figure 3: Flux ratio of the planet over the star as a function of wavelength for 
the two cases considered. Upper panel: no rotation, lower panel: fast rota- 
tion. Several curves using a black-body planet are shown, with temperatures 
representative of the night side and for a planet-averaged stellar heating. 



This planet is very well suited for measuring its light curve 
during a full orbit in order to constrain its heat redistribution 
by having the ther mal curve as a functio n of the orbital phase, 
in a similar way to iKnutson et alJ (120071) . Figure [4] show such 
phase curves for the two scenarios considered. For a non rotat- 
ing planet, as expected the maximum of flux occurs when the 
substellar point is facing the Earth (t — in Figure 0J. On the 



contrary, due to the radiative timescale of the atmosphere when 
a rotation is introduced, the maximum of flux is shifted and we 
measure a delay. In the particular case presented here, since the 
period of rotation is about half the orbital period, the shift is 
about 180°. 

W18b 




Figure 4: Flux ratio of the planet over the star integrated in two Spitzer bandpass 
as a function of time for the two cases considered, t = is when the substellar 
point is facing Earth. Upper panel: no rotation, lower panel: fast rotation. 



4. Conclusions 

We presented a model for the planet WASP- 18b in which the 
stellar heating has been modulated in order to account for the 
different longitude of the planet, with and without rotation. 

In the case with a fast rotation, the temperature mixing is too 
effici ent to match the obse rvational data of secondary eclipse 
from lNvmever et al.l (1201 lh . The day side is too cool indeed to 
get the observed amplitude. 

The no-rotation case, corresponding to a nill heat redistri- 
bution (in the longitud inal direction) fits very well the data of 



iNvmever et al.l (1201 lh . In that case, we expect a full phase 
curve to exhibit a large amplitude and no phase shift with re- 
spect to the substellar point. 

This work does not pretend to investigate all the dynamical 
processes taking place in the exoplanets atmospheres. In partic- 
ular, it neglects any latitudinal heat transport. On the contrary, 
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it aims at showing that a zonal circulation would be in this case 
very low. Moreover, adding latitudinal transport would lead 
to a more uniform global structure. Given the agreement be- 
tween our mod e l with out any transport at all and the data from 



Nvmever et al.l (1201 11). this transport shoul d also be minimal. 



As reviewed in IShowman & Polvanil (1201 lh . in most 3D Gen- 
eral Circulation Models, a zonal superrotating jet at the equator 
is the main feature that develops and can be observed in hot exo- 
planets. In that case our simpler approach can be representative 
of this feature. 

A Spitzer thermal phase cur ve during a full orbit h as already 



been measured for this planet (IMaxted et al.l 12012b . The data 



pertaining to the 3.6 and 4.5 microns light curves agree well 
with the model we present without rotation. Detailed compar- 
ison of the model with the observed light curve will enable us 
to translate the rotation (here seen as sign of heat redistribution) 
into an upper limit of the physical wind speed of the atmosphere 
slow for WASP- 18b. Moreover, an analysis of the shape and 
phase of the peak of the light curve will allow us to constrain 
the heat redeposition processes. 
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